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^ Abstract. We present some mathematical results concerning the strip projection method and a com- 
Cf) puter program for generating quasiperiodic packings of decagonal two shell-clusters. 
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m . 1. The two-shell decagonal cluster C 
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v~j defines the usual two-dimensional representation of the rotation group C5. 
%j By starting from the given vector 



BASIS{:,1) 



BASIS(1,1) 
BASIS(2 1 1) 



we generate the vectors 



BASIS(:, 2) = C5 BASIS(:, 1) 
BASIS{:, 3) = C5 BASIS(:,2) 
BASIS{:,4) = C5 BASIS{:,3) 
BASIS(:,5) = C5 BASIS(:,A). 

The five points of the plane corresponding to 

BASIS(:,1), BASIS{:,2), BASIS{:,3), BASIS{:,4), BASIS(:,5) 
are the vertices of a regular pentagon. If we add the points 

-BASIS{:,1), -BASIS{:,2), -BASIS{:,3), -BASIS(:,A), -BASIS(:,5) 
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we get the vertices of a regular decagon representing the first shell of our two-shell decagonal cluster 
C. 

In a very similar way, by starting from the second given vector 

/ BASIS(1,6) \ 

BASIS(:,6) = 

y BASIS(2, 6) J 

we generate the vectors 

BASIS(:, 7) = C5 BASIS(:, 6) 
BASIS(:, 8) = C5 BASIS(:, 7) 
BASIS(:, 9) = C5 BASIS(:, 8) 
BASIS(:, 10) = C5 BASIS(:, 9). 

The 10 points 

5ASLS(:,6), BASIS(:,7), BASIS{:,8), BASIS{:,9), BASIS(:, 10) 
-5ASLS(:,6), -BASIS(:,7), -BASIS(:,8), -BASIS(:,9), -BASIS(:, 10) 

are the vertices of a regular decagon representing t/ie second shell of our two-shell decagonal cluster C. 



2. The physical two-dimensional space E regarded as a subspace of R 10 

The vectors of M 10 

BASIS(1, :) = ( BASIS(1, 1), BASIS(1, 2), BASIS(1, 3), ... , BASIS(1, 10)) 

BASIS(2,:) = ( 5^/5(2,1), 5^5/5(2,2), BASIS(2, 3), ... , 5ASLS(2, 10)) 
are orthogonal 

10 

( BASIS(1,:), BASIS(2,:) ) = BASIS(1,J) ■ BASIS(2,J) = 

j=i 

and have the same norm 



NORM 
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^ [BASIS(1,J) 
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J2 I BASIS (2, J) 



The subspace E of 1R 10 spanned by these vectors is a two-dimensional space representing the physical 
space. The orthogonal basis 

B = { BASIS(1, :), BASIS(2, :) } 



is very useful when we have to find the orthogonal projection V^V of a vector V G M on E. 
Since 

BASISjl,:) BASIS(2, :) 
NORM ' NORM 



is an orthonormal basis of E we get 

BASIS(1,:) \ BASIS(1,:) I T/ BASIS(2,:) \ BASIS(2,:) 



I BASISQ.,:) \ BASIS(1,:) I BASIS(2, :) \ 
E \ ' NORM J NORM \ ' NORM / 

that is, the coordinates of VeV in the orthonormal basis B are 

y BASISjl,:) \ / y BASIS(2,:] 



NORM 



NORM I ' \ ' NORM 
((V, BASIS(1,:)), (V, BASISp,:))). 



NORM 

In our program we have to find the projection of the points of Z 10 lying in the strip on E. 

In this case the factor 1/NORM is the same for all the points (scaling factor) and we neglect it, that 

is, we use the projector 

V:R 10 — ► K 2 VV = ((V, BASIS(1, :) ), ( V, BASIS(2, :) ) ) . 

4. The window W 

Let V£ : M 10 — ► E- 1 be the orthogonal projector corresponding to the subspace 

E x = { V G M 10 | {V, W} = for all W e E } 

This 8-dimensional subspace (called the internal space) is the orthogonal complement of E. 
The window W is the projection on E -1 

W = 7^(K) 

of the hypercube K = [ -0.5, 0.5 ] 10 , that is, 

K = { ( V(l), 1/(2), ... , V(10) ) | - 0.5 < V{J) < 0.5 for all J G {1, 2, 10} }. 

The window W is a polyhedron in the 8-dimensional space E- 1 . The 7-dimensional faces of W are the 
projections of certain 7-dimensional faces of the hypercube IK from IR 10 . 
Let 

E(I, :) = ( E(1, 1), E(I, 2), E(I, 3), E(1, 10) ) / G {1, 2, 3, 10} 
be the vectors of the canonical basis of M 10 , that is, 

£(1,:) =(1,0,0,0,0,0,0,0,0,0) 
£(2,0 =(0,1,0,0,0,0,0,0,0,0) 
£(3,0 =(0,0,1,0,0,0,0,0,0,0) 



£(10,0 = (0,0,0,0,0,0,0,0,0,1). 



Each 7-face of K is parallel to 7 of these vectors and orthogonal to the other 3 vectors. 
For each three distinct vectors 

£(/i,0, £(A0, £(A0 
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the number of 7-faces of K orthogonal to them is 2 3 . The hypercube IK has 

sets of 2 3 parallel 7-faces. We label them by using the set 

{ (71,72,73) G Z 3 | 1 < 71 < 8, 71 + 1 < 72 < 9, 72 + 1 < 73 < 10 } 
having 210 elements. The set 

[7 G R w 

contains a point and only one from each of the eight 7-faces of IK corresponding to (71, 72, 73). 



77(71), 1/(72), 17(73) G {-0.5, 0.5} 
[7(7) = for 7 G" {71, 72, 73} 



3. The strip § and the quasiperiodic pattern Q 

The strip § is defined as 

§ = { V G M 10 | V\\' G W }. 
In M 3 the cross-product of two vectors v = (v x ,v y ,v z ) and w = (w x .w y ,w z ) can be defined as 



v x w 
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where {i,j, k} is the canonical orthonormal basis of 1R 3 . The vector v x w is a vector orthogonal to v 
and w, and the scalar product between this vector and any vector u = (u x ,u y ,u z ) can be defined as 



(u, v x w) = 



Vx Vy U z 
V X Vy V Z 

W X Wy W Z 



In a very similar way, we obtain a vector W orthogonal to 9 vectors 

U(I, :) = ( C/(7, 1), [7(7, 2), [7(7, 10) ) 7 G {1, 2, 3, 9} 
by expanding the formal determinant 



W = 



75(1,:) 75(2,:) 75(3,:) 

C7(l, 1) C7(l, 2) [7(1,3) 

[7(2,1) [7(2,2) [7(2,3) 

[7(9,1) [7(9,2) [7(9,3) 



75(10,:) 
[7(1,10) 
[7(2,10) 

[7(9,10) 
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containing in the first row the vectors of the canonical basis 5, and 











V(l) 


V(2) 


V(3) 
















5(1,1) 


5(1,2) 


5(1,3) 


... 5(1,10) 








(V,W) = 


5(2,1) 


5(2,2) 


5(2,3) 


... 5(2,10) 














5(9,1) 


5(9,2) 


5(9,3) 


... 5(9,10) 






for any 


vector V from M 10 . 


For example, 
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5(1,1) 


5(1,2) 


5(1,3) 


5(1,4) 


5(1,5) 


5(1,6) 


5(1,7) 5(1. 


8) 5(1,9) 


5(1,10) 




5(2,1) 


5(2,2) 


5(2,3) 


5(2,4) 


5(2,5) 


5(2,6) 


5(2,7) 5(2. 


8) 5(2,9) 


5(2,10) 








5(1,:) 


5(2,:) 


£(3,0 












BASIS(1,1) 


BASIS(1,2) 


5AS5S(1,3) 












BASIS(2, 1) 


BASIS(2,2) 


5AS5S(2,3) 







where 

5(7, J) = BASIS(I, J) 

is a vector orthogonal to the 9 vectors 

5(4,0, £(5,0, £(6,0, £(7,0, £(8,0, £(9,0, £(10,0, BASISQ.,:) and BASIS (2,:). 

Since 5(1, :) - V£E(I, :) is a linear combination of BASIS(1, :) and BASIS(2, :) it follows that W 
is a vector lying in E- 1 orthogonal to the 7 vectors 

7^5(4,0, 7^5(5,0, 7^5(6,0, 7^5(7,:), 7^5(8,0, 7^5(9,:) and 7^5(10, :) 

that is, ortogonal to the two parallel 7-faces of W corresponding to (71, 72, 73) = (1, 2, 3). 
In addition, we get 



(V,W) 



V(l) V(2) V(3) 

BASIS(1,1) BASIS(1,2) BASIS(1,3) 
BASIS(2,1) BASIS(2,2) BASIS(2,3) 



for any 



V = (V(1), V(2), V(3), V(A), V(5), V(6), V(7), V(8), V(9), V(10) ) G 



>io 
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The projection on K 1 - of a point V G Z 10 is "between" the two parallel 7-faces of W corresponding to 
(71, 12, 13) = (1, 2, 3) if and only if 



-STRIP(1, 2,3) < 



V(l) V{2) V{3) 

BASIS(1,1) BASIS(1,2) BASIS(1,3) 
BASIS{2,1) BASIS(2,2) BASIS(2,3) 



where 



STRIP(1,2,3) = 



max 
U(l) e {-0.5, 0.5} 
U{2) e {-0.5, 0.5} 
LT(3) e {-0.5, 0.5} 



< STRIP(1,2,3) 



U(3) 



U(l) U(2) 
BASIS(1,1) BASIS(1,2) BASIS(1,3) 
BASIS(2,1) BASIS(2,2) BASIS(2,3) 



Similar relations can be obtained for each of the 120 systems of indices (11,12,13). 

The quasiperiodic pattern defined by the above construction is 

Q = { VV | V e §nz 10 }. 

We have 

VE{J,:) = {(E{J,:), BASIS(1,:) ), (£(J,:), BASIS(2, :))) = ( BASIS(1, J), BASIS(2, J) ). 

Therefore, for each point VV of Q the nearest and second nearest neighbours are distributed at (some 
of) the vertices of the translation VV + C of the generating two-shell cluster C. The pattern Q can be 
regarded as a quasiperiodic packing of copies of the cluster C. 

More details, bibliography and samples can be found on the website 

http://f pcm5 . f izica.unibuc . ro/~ncotf as 

5. Computer program in FORTRAN and MATHEMATICA 



! QUASIPERIODIC PACKING OF COPIES OF A TWO-SHELL DECAGONAL CLUSTER 

! PLEASE INDICATE HOW MANY POINTS DO YOU WANT TO ANALYSE 
INTEGER, PARAMETER : : N = 10000 

INTEGER I, J, K, L, II, 12, 13, JJ, JP, JPP 
REAL Dl, D2, D3, R, XP, YP 

REAL, DIMENSION(IO) : : V, W, TRANSLATION 
REAL, DIMENSION (2, 10) :: BASIS 
REAL, DIMENSION (2, 2) :: C5 
REAL, DIMENSI0N(1:8,2:9,3: 10) :: STRIP 
REAL, DIMENSION (N ,10) :: POINTS 
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REAL, DIMENSION (N) :: XPOINT, YPOINT 



! PLEASE INDICATE THE COORDINATES OF A POINT BELONGING TO THE FIRST SHELL 
BASIS(l.l) = 1.0 
BASIS(2,1) = 0.0 

! PLEASE INDICATE THE COORDINATES OF A POINT BELONGING TO THE SECOND SHELL 
BASIS(1,6) = 0.9 
BASIS(2,6) = 1.1 

! PLEASE INDICATE THE TRANSLATION OF THE STRIP YOU WANT TO USE 
TRANSLATION =3.7 

C5(l,l) = (SQRT(5.0)-1.0)/4.0 

C5(l,2) = -SQRT(5.0 + SQRT(5 . 0) ) / (2 . * SQRT(2.0)) 
C5(2,l) = SQRT(5.0 + SQRT(5 . 0) ) / (2 . * SQRT(2.0)) 
C5(2,2) = (SQRT(5.0)-1.0)/4.0 
DO J = 2, 5 
DO I = 1, 2 

BASIS(I.J) = C5(I,1) * BASIS(l.J-l) + C5(I,2) * BASIS(2,J-1) 
BASIS(I,5+J) = C5(I,1) * BASIS(1,4+J) + C5(I,2) * BASIS(2,4+J) 
END DO 
END DO 
STRIP=0 
DO II =1, 8 
DO 12 =11+1, 9 
DO 13 =12+1, 10 
DO Dl =-0.5, 0.5 
DO D2 =-0.5, 0.5 
DO D3 =-0.5, 0.5 

R = Dl * BASIS(1,I2) * BASIS(2,I3) + D3 * BASIS(l.Il) * BASIS(2,I2) + & 
D2 * BASISC1.I3) * BASIS (2, II) - D3 * BASIS(1,I2) * BASIS(2,I1) - & 
Dl * BASISC1.I3) * BASIS(2,I2) - D2 * BASIS(l.Il) * BASIS(2,I3) 
IF ( R > STRIP(I1,I2,I3) ) STRIP(I1,I2,I3) = R 
END DO 
END DO 
END DO 

IF( STRIP(I1,I2,I3) . EQ. ) STRIP (II , 12, I3)=N * SUM( BASIS(1,:) ** 2) 
END DO 
END DO 
END DO 

PRINT*, 'COORDINATES OF THE POINTS OF THE TWO-SHELL C5-CLUSTER: ' 
DO J = 1, 10 

PRINT*, J, BASIS(l.J), BASIS(2,J) 
END DO 
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PRINT* , ' * STRIP TRANSLATED BY THE VECTOR WITH COORDINATES : ' 
PRINT* , TRANSLATION 

PRINT*, '* PLEASE WAIT A FEW MINUTES OR MORE, & 

DEPENDING ON THE NUMBER OF ANALYSED POINTS' 

JP = 

POINTS = 
POINTS (1,:) = ANINT( TRANSLATION) 
K = 1 
L = 

DO I = 1, N 

V = POINTS (I, : ) - TRANSLATION 
JJ = 1 
JPP = 

DO II =1, 8 
DO 12 =11+1, 9 
DO 13 =12+1, 10 

R = V(I1) * BASIS(1,I2) * BASIS(2,I3) + V(I3) * BASIS(1,I1) * BASIS(2,I2) + & 
V(I2) * BASIS(1,I3) * BASIS(2,I1) - V(I3) * BASIS(1,I2) * BASIS(2,I1) - & 
V(I1) * BASIS(1,I3) * BASIS(2,I2) - V(I2) * BASIS(1,I1) * BASIS(2,I3) 

IF ( R < - STRIP(I1,I2,I3) .OR. R > STRIP (II , 12 , 13) ) JJ = 
IF ( R == - STRIP(I1,I2,I3) .OR. R == STRIP(I1 , 12 , 13) ) JPP = 1 
END DO 
END DO 
END DO 

IF( J J .EQ. 1 ) THEN 
XP = SUM( V * BASISC1, :) ) 
YP = SUM( V * BASIS (2, : ) ) 
13 = 1 
DO J = 1, L 

IF( XP == XPOINT(J) .AND. YP == YPOINT(J) ) 13 = 
END DO 

IF( 13 == 1) THEN 

IF( JPP .EQ. 1 ) JP = JP + 1 

L = L + 1 

XPOINT(L) = XP 

YPOINT(L) = YP 

ELSE 
END IF 

DO II = 1, 10 
DO 12 = -1, 1 

W = POINTS (I, : ) 
W(I1) = W(I1) + 12 
13 = 
DO J = 1, K 

IF( ALL(W .EQ. POINTSCJ,:)) ) 13 = 1 
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END DO 

IF ( 13 == .AND. K < N ) THEN 

K = K + 1 

POINTS (K, : ) = W 

ELSE 

END IF 
END DO 
END DO 
ELSE 
END IF 
END DO 

PRINT*, ' NUMBER OF ANALYSED POINTS :', K 

PRINT*, 'NUMBER OF OBTAINED POINTS : ' , L 
PRINT*, 'NUMBER OF POINTS LYING ON THE FRONTIER OF THE STRIP:', JP 
PRINT*, 'PLEASE INDICATE THE NAME OF A FILE FOR RESULTS' 
WRITE(4,10) 

10 FORMAT ( ' Show [Graphics [{PointSize [0 . 03] , { ') 
DO J = 1, 10 

WRITE(4,20) BASIS(l.J), BASIS(2,J) 

20 FORMAT ( ' Point [{ 'F10 . 5 ' , ' F10 . 5 , ' }] , ') 

END DO 

DO J = 1, 10 

WRITE(4,30) -BASIS(1,J), -BASIS(2,J) 

30 FORMAT ( ' Point [{ 'F10 . 5 ' , ' F10 . 5 , ' }] , ') 

END DO 

WRITE(4,40) 0.0, 0.0 

40 FORMAT ( ' Point [{ 'F10 . 5 ' , ' F10 . 5 , ' }] } }] , & 

PlotRange -> All, AspectRatio -> 1 ]') 

WRITE (4, 80) 

80 FORMAT ('Show [Graphics [{PointSize [0 . 02] , { ') 
DO J = 1, L-l 

WRITE(4,90) XPOINT(J), YPOINT(J) 

90 FORMAT ( ' Point [{ 'F10 . 5 ' , ' F10 . 5 , ' }] , ') 

END DO 

WRITE(4,100) XPOINT(L), YPOINT(L) 
100 FORMAT ( 'Point[{'F10.5' , 'F10.5, '}]} } ], & 

PlotRange -> All, AspectRatio -> 1]') 
PRINT*, '* OPEN THE FILE CONTAINING THE RESULTS WITH "NotePad" ' 
PRINT*, '* SELECT THE CONTENT OF THE FILE ("Select All") & 

AND COPY IT ("Copy")' 
PRINT*, '* OPEN "MATHEMATICA", PASTE THE COPIED FILE, & 

AND EXECUTE IT ( "Shif t+Enter" ) . ' 

END 
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